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The existence of unique scaling in a crossover regime between viscous and inertial hydrodynamic 
regimes is revealed for homogeneous, isotropic, incompressible, spinodal turbulence which is char- 
acterized, to begin with, by three different length scales and a velocity scale. The obtained scaling 
exponents are found to be in agreement and in consistency with available simulation results for a 
broad range of crossover regime. Also, it is observed that the spinodal turbulence in the crossover 
regime is in complete consistency with the universality class of self-preservation of decaying grid 
turbulence. We then obtain analytical forms for various scalings, valid in the crossover regime, 
through the analysis for self-preservation of spinodal turbulence. 
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A binary fluid homogeneous mixture, when quenched 
to thermodynamically unstable state, starts separating 
with growing size of domains of two different components 
of the binary fluid. Here our focus is on incompressible 
50/50 mixture of two fluids having identical viscosity, 
density and undergoing the separation process. Through 
experimental, theoretical and simulation work, a general 
understanding has evolved on various scaling phenom- 
ena which unfold during the course of the separation 
or spinodal decomposition process (see refs. Q, Hi HI 
and references cited therein). These scaling phenomena 
can be categorized into four different regimes, namely, 
1) pure diffusive, 2) viscous hydrodynamic, 3) inertial 
hydrodynamic, and 4) crossover regime between viscous 
and inertial hydrodynamic. Theoretically predicted scal- 
ings L{t) t^/"^ and L{t) t of growing domain size L 
in time t during the early stage of separation in diffusive 
and viscous hydrodynamic regimes j^], respectively, are 
well supported (see refs. 0, 0, ^^'^ references cited 
therein). The pioneer scaling L{t) ^ t^^^ in the later 
stage of inertial hydrodynamic regime 0, though has 
found support 0, still remains controversial due to 
the physical requirement of Reynolds number saturation 
in the long time limit and which has led to exponent val- 
ues other than 2/3 [i0,E3|- 

Contrary to this, we also 
note that scaling theory of Kendon flil obtains 2/3 and 
is consistent with Reynolds number saturation. Further, 
a broad range of crossover between viscous and inertial 
hydrodynamic regimes exists as exhibited in simulation 
studies 1^ 1^ . It is clearly apparent from the results of 
these simulation studies that scaling does exists in the 
crossover regime, but so far it has not been dealt with 
and predicted theoretically. In this letter, we reveal ex- 
istence of unique scaling L ^ t'^/'* in the broad crossover 
regime, connecting viscous and inertial hydrodynamic 
regimes where viscous and inertial effects are both im- 
portant. This scaling along with other obtained scalings 
for velocity and its length scales, Taylor-scale and Kol- 



mogorov dissipation scale are found to be in agreement 
with the simulation data of Kendon et al. for spin- 
odal turbulence and in complete consistency with results 
for decaying of homogeneous, isotropic turbulence jl3l |. 
as well. For our analysis, we use theoretical framework 
of Kendon 1^ . 

The Navier-Stokes equation governing the phenomena 
of spinodal turbulence velocity field v in incompressible, 
isothermal 50/50 mixture of two fluids of density p can 
be written as 



dv 

dt 



p(v • V)v = -Vp + ryV^v - Vr. 



(1) 



Here p is pressure, rj is viscosity and —VP is the inter- 
facial force density, usually approximated hy a/L^ 0,0 
where a is the interfacial tension between the two flu- 
ids and L is length scale representing the domain size. 
For this homogeneous isotropic case, the kinetic energy 
of turbulence per unit volume, i.e. ^(v • v), is governed 
by 



dt 



p{v ■ v)/2 = -r?((Vv)2} + (-V ■ VP) 



(2) 



where ( ) denotes ensemble average, the first term on the 
right hand side (rhs) represents viscous dissipation rate 
and the second source term on the rhs accounts for the 
energy transfer rate to the fluid arising from the interfa- 
cial force. This source term is approximated, so far, by 
aL/L'^ 12] for the scaling purpose and implicitly assumes 
the approximation w ^ L ^ L/t. We should mention 
that we do not use such an approximation for v when 
a separate scaling (see scaling given by Q) is used for 
it. In this respect we allow Navier-Stokes equation, gov- 
erning the phenomena, to govern and yield length scale 
Ly ^ \v\t (or V ^ Ly) of the velocity. 

Now consider the scaling behavior of various quanti- 
ties, namely, three length scales L, , L^2 and velocity 
scale V, as 



domain size: 



(3) 
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fluid velocity: y ^ t*^ , 
velocity first derivative: Vv ~ - — ^ i*^ " 



velocity second derivative: V v 



7^2 



(4) 
(5) 



Further, using scalings L ^ t", v ^ and approxima- 
tion —VP ~ cr/i^, we can write for the source term 



(-V • VV) - a 



L2 



at 



(7) 



and which is different than the usually employed approx- 



imation aL/L ^ at 



by Kendo 



f empi( 

n Ea. 



We should 



mention that this new relation Q leads us to a unique 
set of values for all scaling exponents (a, /3, a', a") and 
yields a situation where scaling exponents of all terms 
in the Navier-Stokes equation ^ are identical. For this 
reason, we have associated the obtained unique scaling 
to the crossover regime, between viscous and inertial hy- 
drodynamic regimes, in which all terms in the Navier- 
Stokes equation are relevant to the involved scaling phe- 
nomenon. 

Using scaling behaviors given by l|SIl-|3l allow us to 
write the turbulent kinetic energy Eq. ([SJ in terms of 
power of t as 



+ at^'- 



and which further suggests 



This gives 



2/3 - 1 = 2/3 - 2a' = /3 - 2a. 



1/2, P=l-2a. 



(8) 



(9) 



(10) 



Now we study the behavior of Navier-Stokes equation 
when various scalings ®-® are substituted in it. It 
should be noted that in the absence of —VP, scaling of 
— Vp is identical to the scaling of convection term, i.e., 
second term on the left hand side of Eq. . We make use 
of this fact along with and approximation —VP ^ 

a/L"^ to obtain, from the Navier-Stokes equation, 



Further making use of (|10|) into (|ll|l . we obtain 



pl3t 



-2a 



pr^ 



which suggests 



- 2a = -4a + 3/2 = 1 - 2a - 2a" 
and yields a unique solution 

a = 3/4, a" = 1/2. 



(11) 



(12) 



(13) 



(14) 



Subsequently from (|10|l we obtain 
/3--1/2. 



(15) 



This unique scaling solution suggests that all the terms, 
including pressure term, in the Navier-Stokes equation 
are scaled identically as t~^/'^. For the purpose of com- 
pleteness, we must mention that the obtained unique set 
of scaling exponents (a = 3/4,/? — —1/2, a' = 1/2, a" = 
1/2) along with — v • VP ^ aj^ at^~'^°' also provide 
identical scaling exponent for all the terms in the equa- 
tion for local energy or instantaneous kinetic energy 
of turbulence pv ■ v/2. This further strengthen the as- 
sociation of obtained scalings to the crossover regime in 
which all the terms in the equation for local energy be- 
come crucially important. 

A few more scaling relations for velocity length scale 
Ly, two different Reynolds numbers {Rcl, Rsv), dissi- 
pation rate e, Taylor-scale A and Kolmogorov dissipation 
scale Xd can be derived using the above scaling exponents. 
These are 



ReL 

i?ev 

e 

A 

Ad 



,1/4 



|v|t^fl/2 

Pl(v V)v| ,0 
7y|V2v| ^ ' 



,1/2 



„3 \ 1/4 

V ^ ^1/2 



(16) 
(17) 

(18) 
(19) 
(20) 

(21) 



It should be noted that scaling for velocity length scale 
Ly ^ t^/^ turns out to be different than the scaling for 
domain size L ^ t^/^. In fact, this prediction is consistent 
with simulation study by Kendon et al. as exhibited 
in their figure (11a) in the crossover regime 10^ < t < 
10^. Thus in view of this, the traditional approximation 
aL/L^ for (— v • VP) which employs |v| ~ L is not valid 
in the crossover regime. And our choice of approximation 
as given by turns out to be correct for the crossover 
regime. 

Now important question about the validness of var- 
ious obtained scalings arises. The available simulation 
study (hereafter referred to as KCPDB paper) pro- 
vides information for a long range of crossover regime 
(10^ < t < 10^ or 1 < ReL < 50 ). In KCPDB, results 
for L as shown in figure 6 and values of a listed in Ta- 
ble 5 suggest that our predicted exponent a = 3/4 is a 
very reasonable value in the crossover region. Result of 
figure (17a) in KCPDB provide scaling close to for 
dissipation rate and is in agreement with present scaling 
given by H19|) . It is encouraging to note that the pre- 
dicted scaling exponent 1/2, given above in (|20|) and H21|l . 
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for Taylor-scale and Kolmogorov dissipation scale are in 
agreement with crossover regime data of figure (18b) in 
KCPDB. Further, the present scaling prediction t"-^/^ for 
each term of Navier-Stokes equation appears visually to 
be consistent with data of figure (18a) in KCPDB. 

Apart from the consistency with spinodal turbulence 
simulation data, we note that the obtained scalings are 
also consistent with decaying, homogeneous isotropic tur- 
bulence. In this case, when the turbulence velocity scale 
is decaying sls v ^ with /3 as a constant, A ~ t^/^ 
and integral length scale of velocity t^^^ (see e.g. refs. 
[l^ll3 |\ So these exponents for A and the integral length 
scale of velocity are exactly the same as predicted above 
if we associate L„ to the integral length scale. In addition 
to this, the above predicted scaling exponent for veloc- 
ity decay (3 = —1/2, scalings of dissipation rate ^ 
and Kolmogorov dissipation scale ~ t^/^ are identical to 
those belonging to universality class of initial period of 
decay of grid turbulence preserving the shape of whole of 
the energy spectrum and correlations functions, except 
for the portion at very smallest values for magnitude of 
wave numbers |l3l |. In view of this consistency with the 
universality class, we can conclude that spinodal turbu- 
lence is self-preserving in the cross over regime. Now 
we provide analysis on self-preservation of spinodal tur- 
bulence which obtains exact analytical forms for various 
scalings consistent with the scaling exponents obtained 
above. 

We should note that though we started with three 
length scales, the scaling exponents for Lv and L^2 are 
identical. This suggests reduction by one in number of 
length scales to characterize crossover regime of spinodal 
turbulence. For the analysis of self-preservation, we con- 
sider equation for 



i?,,(r, t) = (v(x, t) ■ v(x + r, t)) = 3fir, t) + r 



dr 
(22) 

for homogeneous iso trop ic turbulence which can be ob- 
tained from Eq. [lj| and can be written as 

9R^^ d \ f d 4\ , , , 

/ ^2 2 d\ 

7r^ + -7r rR" + 2(v(x,<).s(x + r,t)). (23) 
\ or'^ r or / 

Here v? = (v • v)/3, Rij is two-point velocity correla- 
tion, r — |r| is distance between two points for two- 
point correlation, / is the longitudinal velocity corre- 
lation coefficient, k{r,t) is single scalar function deter- 
mining the triple- velocity correlations (see ref. ^3|) and 
s = —VP represents interfacial force per unit volume. 
Also, for the homogeneous isotropic turbulence we can 
write (v(x, i) •s(x + r, i)) = usT{r, t) where = (s •s)/3. 

Now we consider characteristic length scale I = l{t) 
(which will turn out to be similar to Lv), introduce a 



dimensionless variable ip = r/l and take self-preserving 
form for functions /, fc, as 



Substituting these in Eq. H23I) . we obtain 

^^^^'^^'^^di^^'^^ " ^^3(^) + 

-I- usJ-(ip) 



(24) 



(25) 



where Fi, F2, F3 and F4 are some functions of ip only and 
their details are not presented here as they are not rele- 
vant for further analysis. For self-preservation, Eq. H25|l 
suggests that all ratios between different factors multi- 
plying Fi, i^2, ^3, F4 and T should be constants. For ob- 
taining analytical solutions, consider 



pu^ dl 
~rdi 

~T 

du^ 



A 



= A- 



rju 



, pu^ dl 



dt 



us 



A, 



pu 

T 



I dt 
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(26) 
(27) 
(28) 
(29) 



where Ai,A2, A3 and A4 are constants. These Eqs. (|26|) - 
(|29|l suggest analytical solutions, written as 



2Al7? 



{t - to) + ll 



1/2 



A^Ji 
P 



2Ai7] 



{t - to) + ll 



^1/2 



A3 



-1/2, 



s^^{A,^f 
P 



2Air/ 



-1 -3/2 



[t - to) + ll 



(30) 



(31) 



(32) 



(33) 



Here Iq is value of length scale I at time to which is in the 
crossover regime. It should be noted that I ~ t^^"^ and 
thus is similar to Lv- Also, s ^ t^^/^ is consistent with 
the scaling L ^ t^/^ as s '~ a /L"^ . 

In conclusion, we have suggested existence and values 
of the unique scaling exponents in the crossover regime 
for the phenomena of spinodal turbulence where all forces 
in the Navier-Stokes equation are important and scale 
to identical exponent. The obtained scalings have been 
found to be in reasonable agreement with the simula- 
tion data and is in complete consistency with the self- 
preserving phenomena of decay of grid turbulence dur- 
ing the initial period. Further, analytical forms for char- 
acteristic length, velocity scales (/ and u, respectively) 



and rms of interfacial force density component, i.e. s, 
are obtained through the analysis for self-preservation 
of spinodal turbulence in the crossover regime. Though 
there exists a few studies on the crossover regime 0, ^| , 
it is hoped that the present work would stimulate fur- 
ther more focused studies on the scaling behaviour of 
crossover regime of the spinodal turbulence so as to rig- 
orously verify or dispute the scaling exponents predicted 
theoretically in this letter. 

I am very much thankful to my friend Dr. Paul Stansell 
for introducing me to this subject of binary fluid demix- 
ing. 
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